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ABSTRACT 

Radiation pressure acting on gas and dust causes H II regions to have central den- 
sities that are lower than the density near the ionized boundary. H II regions in static 
equilibrium comprise a family of similarity solutions with 3 parameters: /3, 7, and the 
product Qo^-rms; P characterizes the stellar spectrum, 7 characterizes the dust/gas ra- 
tio, Qo is the stellar ionizing output (photons/s), and n rms is the rms density within 
the ionized region. Adopting standard values for /3 and 7, varying Qo^rms generates 
a one-parameter family of density profiles, ranging from nearly uniform-density (small 
<5o n rms) to shell-like (large Qon rms ). When Qon rrns ^ 10 52 cm~ 3 s _1 , dusty H II regions 
have conspicuous central cavities, even if no stellar wind is present. For given f3, 7 and 
Qo n rms> a fourth quantity, which can be Qo, determines the overall size and density of 
the H II region. Examples of density and emissivity profiles are given. We show how 
quantities of interest - such as the peak-to-central emission measure ratio, the rms- 
to-mean density ratio, the edge-to-rms density ratio, and the fraction of the ionizing 
photons absorbed by the gas - depend on /3, 7, and Qonrms- For dusty H II regions, 
compression of the gas and dust into an ionized shell results in a substantial increase in 
the fraction of the stellar photons that actually ionize H (relative to a uniform density 
H II region with the same dust/gas ratio and density n = n rms ). We discuss the extent 
to which radial drift of dust grains in H II regions can alter the dust-to-gas ratio. The 
applicability of these solutions to real H II regions is discussed. 

Subject headings: ISM: bubbles; dust, extinction; HII regions; ISM: structure; infrared: 
ISM; radio continuum: ISM 



1. Introduction 



Stromgren! (j 19391 ) idealized photoionized nebulae around hot stars as static, spherical regions 



with a uniform density of ionized gas out to a bounding radius R. The Stromgren sphere model 
continues to serve as the starting point for studies of H II regions around hot stars. However, a 
number of physical effects lead to departures from the simple Stromgren sphere model: dynamical 
expansion of the H II region if the pressure in the surrounding neutral medium cannot confine the 
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ionized gas; deviations from sphericity due to nonuniform density; motion of the star relative to the 
gas; injection of energy and momentum by a stellar wind; absorption of H-ionizing photons by dust 
grains; and radiation pressure acting on gas and dust. Each of these effects has be e n the object of 
a number of investigations, beginning with the study of ionization fronts by iKahnl (| 19541 ). 



Savedoff k, Greend ()1955l ) appear to have been t he first to discus s the expans ion of a spherical 



H II r egion in an initially uniform neutral medium. iMathewd (|l967l . 1 19691 ) and iGail &: Sedlmavr 
(|1979l ) calculated the dynamical expansion of an H II region produced by an O star i n a medium that 
was i nitially neutral, including the effects of radiation pressure acting on the dust. iMathewsl (|1967l . 
19691 ) showed that radi ation pressure on dust woul d produce low-density central cavities in H II 
regions. More recently, iKrumholz &: Matznerl (|2009i ) reexamined the role of radiation pressure on 
the expansion dynamics of H II regions, concluding that radiation pressure is generally unimportant 
for H II regions ionized by a small number of stars, but is important for the expansion dynamics of 
giant H II regions surrounding clusters containing many O-type stars. Their study concentrated on 
the forces acting on the dense shell of neutral gas and dust bounding the H II region, hence they 
did not consider the density structure within the ionized region. 

Dust absorbs hu > 13.6 eV photons that would oth e rwise be able to ionize hydrogen, thereby 
reducing the extent of the ionized zone. iPetrosian et al.l (119721 ) developed analytic approximations 
for dusty H II regions. They assumed the gas density to be uniform, with a constant dust-to-gas 
ratio, and found that dust could absorb a substantial fraction of the ionizing photons in dense H II 
regions. Petrosian et al. did not consider the effects of radiation pressure. 



Dopita et al.l (|2003l . 120061 ) constructed models of compact H II regions, including the effects of 
radiation pressure on dust, and presented models for different ionizing stars and bounding pressures. 
In these models, radiation pressure produces a density gradient within the ionized gas. 

The present paper provides a systematic discussion of the structure of dusty H II regions 
that are assumed to be in equilibrium with an external bounding pressure. The assumptions and 
governing equations are presented in $2j where it is shown that dusty H II regions are essentially 
described by a 3-parameter family of similarity solutions. In <j3]we show density profiles for selected 
cases, as well as surface brightness profiles. The characteristic ionization parameter Ui/ 2 and the 
fraction (1 — /j on ) of the ionizing photons absorbed by dust are calculated. Dust grain drift is 
examined in where it is shown that it can alter the dust-to-gas ratio in the centers of high 
density H II regions. The results are discussed in §5\ and summarized in $6j 



2. Equilibrium Model 

Consider the idealized problem of a static, spherically-symmetric equilibrium H II region ion- 
ized by a point source, representing either a single star or a compact stellar cluster. Assume a 
constant dust-to-gas ratio (the validity of this assumption will be examined later). For simplicity, 
ignore scattering, and assume a^, the dust absorption cross section per H nucleon, to be independent 
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of photon energy hv over the ~ 5eV to ~ 30 eV range containing most of the stellar power. 

Let the star have luminosity L = L n + L. L = L39IO 39 ergs s , where L n and Lj are the luminosi- 
ties in hv < 13.6 eV and hv > 13.6 eV photons, respectively. The rate of emission of hv > 13.6 eV 
photons is Qo = 10 49 <3o,49 s _1 and the mean energy of the ionizing photons is (hv) ; = Lj/Qp. A 



singl e main sequence star of spectral type 06V has L39 = 0.80 and Qo,49 = 0.98 (jMartins et al 



2005). A compact cluster of OB stars might be treated as a point source with much larger values 
of Qo,49 and L39. 

Ignore He, and assume the H to be nearly fully ionized, with photoionization balancing "Case 
B" radiative recombination, with "on-the-spot" absorption of hv > 13.6 eV recombination radia- 
tion. Take the effective radiative recombination coefficient to be as ~ 2.56 x 10 _13 T^~ 0-83 cm 3 s _1 
for 0.5 < T4 < 2, with T 4 = T/10 4 K, where T is the gas temperature. 

Assume the gas to be in dynamical equilibrium (the neutral gas outside the ionized zone is 
assumed to provide a confining pressure). Static equilibrium then requires that the force per unit 
volume from radiation pressure be balanced by the pressure gradient: 

[L n e- T + Li(j)(r)\ 2 ( hu )i d , n , ^ n m 

na d — k h Qbi r (2nfcT) = , (1) 

47rr z c c dr 

where n(r) is the proton density, Lj0(r) is the power in hv > 13.6 eV photons crossing a sphere 
of radius r, and r(r) is the dust absorption optical depth. Eq. (fTJ) underestimates the radiation 
pressure force, because it assumes that recombination radiation (including Lyman-a) and cooling 
radiation escape freely. 

The functions <f>(r) and r(r) are determined by 

d(f) 1 o n 

— = -— asn 47rr - na d <p , (2) 
dr Q 

-j- = n &d j (3) 
dr 

with boundary conditions (f)(0) = 1 and t(0) = 0. Define a characteristic density and length scale 

,16 <?0,49 f (J}V_ 

T?- 83 Vl8eV 
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Aq ^ _^(r£^l) = 2 .47xl0 1B ^ l^fcm , (5) 



and the dimensionless parameters 



„ , ^.^-LMT^-f"*)-! , (6) 
Li Li Q ,49 V \hv)i J 
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The parameter (3, the ratio of the power in non-ionizing photons to the power in photons 
with hu > 13.6 eV, depends solely on the stellar spectrum. We take j3 = 3 as our standard 
value, corresponding to the spectrum of a T* = 32000 K blackbody, but we also consider (3 = 2 
(T* = 45000 K) and (3 = 5 (T* = 28000 K); the latter value might apply to a cluster of O and B 
stars. 
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Fig. 1. — (a) Absorption cross section per H, and (b) radiation pressure cross section per H, averaged over bl ackbody 



igcd c 

spe ctra, as funct i ons of the blackbody temperature for the dust models of IWeingartner fc Draind (|200ll . WD01) 



and lZubko et al.l (|200J, ZDA04). Broken lines show averages over hv < 13.6 eV only, appropriate for dust in neutral 
gas. 



Momentum can be transferred to a dust grain by photon absorption, but also by scattering. 
The cross section ad appearing in eq. ([I]) should be (a pT ), the radiation pressure cross section per 
H, cr pr (^) = cr a bs + (1 — (cos 9) )ci sca , averaged over the spectrum of the radiation field at radius r, 
where cr a b s (^) and a sca (v) are the absorption and scattering cross section per H, and (cos 9) is the 
mean value of the cosine of the scattering angle 9 for photons of frequency v. 

In eqs. ([2]) and ([3|), ad characterizes the effectiveness of the dust in attenuating the radiation 
field. While scattering does not destroy the photon, it does increase the probability of the photon 
undergoing subsequent absorption. Thus, the value of ad in eqs. ([2]) and ([3]) should exceed (<7 a b s )- 

Figure [U(a) shows the dust absorption cross section per H nucleon averaged over a blackbody 
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spectrum, for two dust models (jWeingartner fc Drainei l200ll ; IZubko et al.l 120041 ) that reproduce 
the wavelength-dependent extinction in the diffuse interstellar medium using mixtures of PAHs, 
graphite, and amorphous silicate grains. Fig. [I^b) shows that (<7 pr ), the radiation pressure cross 
section averaged over blackbody spectra, is only slightly larger than (cr a bs)- Given the uncertainties 
in the nature of the dust in H II regions, it is reasonable to ignore the distinction between (cr pr ) 
and the attenuation cross section and simply take ad = (<7 pr ) in eq. (pj|3|). 

For dust characteristic of the diffuse ISM, one could take (<r pr ) ~ 1.5 x 10~ 21 cm 2 FT 1 for 
2.5 x 10 4 K < T rac j < 5 x 10 4 K. However, dust within an H II region may differ from average 
interstellar dust. For example, the small-size end of the size distribution might be suppressed, in 
which case would be reduced. Low metallicity galaxies will also have lower values of ad, simply 
because there is less material out of which to form grains. In the present work we will assume a 
factor ~1.5 reduction in ad relative to the local diffuse ISM, taking ad ~ 1 x 10~ 21 cm 2 H^ 1 as the 
nominal value, but larger and smaller values of ad will also be considered. 

The dimensionless parameter 7 (defined in eq. [7j) depends also on the gas temperature T and 
on the mean ionizing photon energy {hv)i, but these are not likely to vary much for H II regions 
around OB stars. We take 7 = 10 as a standard value, but will also consider 7 = 5 and 7 = 20. 
Low-metallicity systems would be characterized by small values of 7. 

Switching to dimensionless variables y = r/Ao, u = no/re, the governing equations {IMS]) become 



du 
dy 
d(j) 
dy 
dr 
dy 



1 - 7 {Pe- T 

y 2 (j) 

— -7- , 

it z u 



7 
u 



u 
T7i 



(8) 
(9) 
(10) 



with initial conditions (p(0) = 1 and r(0) = 0. The solutions are defined for < y < y maX ) where 
2/max is determined by the boundary condition </>(y max ) = 0. The actual radius of the ionized zone 
is R = ymaxAo- For each solution u{y) the mean density is 



3n r*™ 1 2 
W = — — / -y dy , 

2/max JO u 



(11) 



the root-mean-square density is 



^rms — ^0 



3 ^max 1 ■ 

-3 — / — V dy 

2/max JO u 



1/2 



(12) 



and the gas pressure at the edge of the H II region is 



Pedge = 2n(R)kT 



2n kT 

u(y max ) 



(13) 
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Let 

R s o - ( -r^-) = 2.10 x 10 18 ^f T 4 °" 28 cm (14) 
V 47rn rms as / n ^ g 

be the radius of a dustless Stromgren sphere with density n rms = 10 3 n rmSj 3 cm -3 . The fraction of 
the hv > 13.6 eV photons that are absorbed by H is simply 

, Rso ) ^ ^ 

For given (/3,7), varying the initial valued of u = no/n at some fixed y = t/Xq generates 
solutions with different density profiles. Therefore the full set of solutions forms a three-parameter 
family of similarity solutions u(y), cp(y), and r(y), parametrized by /?, 7, and a third parameter. 
The third parameter can be taken to be Qon Tms . For dusty H II regions, an alternative choice for 
the third parameter is the dust optical depth on a path R s o with density Ti rms : 

Td,o = n rms R s0 a d = 2.10 (Q ,49n rmSi3 ) 1/3 T^ 28 ——^ — ? (16) 

10 Z1 cm z 

= 0.188 7 (Q 0i49 „ rinSi 3)V3 rr i.55 i^i (17) 



The static H II regions described by eq. (p]-[3|) are determined by 7 distinct dimensional 
quantities: three parameters describing the central star (Qo, (hu)i, and L n ), the recombination 
rate coefficient ob, the thermal energy kT, the dust cross section per nucleon ad, and the external 
pressure p e dgc confining the H II region. According to the present analysis, this 7-parameter family 
of solutions actually reduces to a 3-parameter family of similarity solutions. The dimensionless 
parameters f3 and 7, plus choice of an initial value for the function u near y = 0, suffice to 
determine the scaled density profile n(r)/no and radius y max = R/^o : this is the 3-parameter 
family of similarity solutions. 

Specifying numerical values for the ratios Qo/ob and kT/(otB{hv)i) fixes the values of uq and 
Ao, thus giving n(r) for r < R. Thus far we have invoked 5 independent parameters, but have not 
actually specified either kT or a^. 

Specifying kT and a_e - the 6th and 7th parameters - allows us to compute the actual values 
of Qo and {hv)i, and the bounding pressure p e dge = 2n(R)kT. 

If the "initial value" of u near the origin is taken as a boundary condition, then p e dge emerges 
as a derived quantity. However, if we instead regard p e dge as a boundary condition, then the initial 
value of u ceases to be a free parameter, and instead must be found (e.g., using a shooting technique) 
so as to give the correct boundary pressure p e dge : the initial value of u is thus determined by the 
7 physical parameters (Qo, {hv)i, L n , as, T, ad, and p e dge)- 



For 7 > 0, u oc exp[(/3 + 1)7/3/] ~* 00 as V ~ ► 0, and the integration must start at some small y > 0. 
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Thus we see how the 3 parameter family of dimensionless similarity solutions corresponds to 
a 7 parameter family of physical solutions. 




Fig. 2. — Normalized density profiles of static equilibrium H II regions, as a function of r/R, where R is the radius 
of the ionized region. Profiles are shown for 7 values of Qo^rms; the numerical values given in the legends assume 
T 4 = 0.94 and (hv)i = 18 eV. (a) Dustless (7 = 0); (b) 7 = 5; (c) 7 = 10; (d) 7 = 20. 



3. Results 

Equations ([8l fl"0j) can be integrated numerically. Figure [2^, shows the solution for the case 
where no dust is present (7 = 0). Radiation pressure associated with photoionization produces a 
density gradient in the H II region, but it is modest unless Qo^rms is very large. The central density 
is nonzero. For Qo,49 n rms 10 3 cm -3 , the density is uniform to within ±15%. 

As discussed above, the dust abundance relative to H is characterized by the parameter 7. 
Density profiles are shown in Fig. [2b-d for f3 = 3 and 7 = 5, 10, 20, corresponding approximately to 
ad = 0.5, 1, 2 x 10~ 21 cm 2 . With dust present, the density formally goes to zero at r = 0. For fixed 
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Fig. 3. — Ionization parameter U1/2 at the half-ionization radius in dusty H II regions (see text) , calculated assuming 
T 4 = 0.94 and (hv)i = 18 eV. 



7, the size of the low-density central cavity (as a fraction of the radius R of the ionization front) 
increases with increasing Qo n rms- The enhancement of the density near the ionization front also 
becomes more pronounced as Qo^rms is increased. For f3 = 3, 7 = 10 and <5o,49"-rms = 10 5 cm~ 3 , 
we find n(R) = 2.5n rms . 

The state of ionization of the gas is determined by the hardness of the radiation field, and the 
value of the dimensionless "ionization parameter" 

U EE H{hU > /H) , (18) 

where n(hu > Ih) is the density of photons with hv > Jh- Within an H II region, the value of 
U varies radially. As a representative value, we evaluate U1/2, the value at the "half-ionization" 
radius R1/2, the radius within which 50% of the H ionizations and recombinations take placed] In 



2 Some authors (e.g., iDopita et al. 200a) use the volume-averaged ionization parameter (U)v- For a uniform 
density dustless HII region, (U)v = (81/2567r) 1/3 (Qs /3 /c)(Qon rms ) 1/3 = 2.83 U 1/2 - 
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(a) Dustless H II Regions 

Emission Measure Profile 




(c) Dusty HII Regions with = 3, 7=10 
Emission Measure Profile 
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(b) Dusty HII Regions with jS=3, 7=5 
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(d) Dusty HII Regions with /3 = 3, 7 = 20 
Emission Measure Profile 
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Fig. 4. — Normalized emission measure (EM) profiles for a cut across the center of H II regions with (a) 7 = (no 
dust), (b) 7 = 5, (c) 7 = 10, and (d) 7 = 20. Profiles are shown for selected values of Qo^rms- Numerical values of 
Qo,49«-rms assume T4 = 0.94 and (hv)i = 18 eV. 



a uniform density dustless H II region, Ri/ 2 = 2 1 / 3 i?5o is the same as the half-mass radius, and 

2/3 

r/ (nodust) _ a B (q a/3 /-.qn 

U l/2 ~ ( 72vr )l/3 c WO^rmsJ • (W) 

For our present models, 

^1/2 = "j-^2 2 ' V 20 ) 

4vrA5n c y{ /2 

where y\i% = -R1/2M0 is the value of y within which 50% of the H ionizations and recombinations 
take place. Figure [3] shows £/i/ 2 as a function of Qon rras for static dusty H II regions with radiation 
pressure, for selected values of f3 and 7. For small values of Qo n rms, dust and radiation pressure 
are negligible and [/i/ 2 coincides with Uy^^ (eq. IT9]) . However, for large values of Qon rms , Uy 2 

falls below ■ F° r 7 ~ 10 - corresponding to the dust abundances that we consider to be 
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likely for Galactic H II regions - we see that U1/2 ~ 0.07 ± 0.02 for Qo,49"rms ^ 10 4 cm 3 . 

The emission measure EM(b) = J v? e ds is shown as a function of impact parameter b in 
Figure [H For small values of Qo n rms) the intensity profile is close to the semicircular profile of a 
uniform density sphere. As Qon rms is increased, the profile becomes flattened, but, if no dust is 
present (7 = 0, Fig. [2^i) , the ionized gas only begins to develop an appreciable central minimum 
for Qo,49n rm s > 10 4 - 5 cm -3 . 

When dust is present, however, the profiles are strongly affected. For standard parameters (3 = 
3, 7 = 10, the emission measure shows a pronounced central minimum for Qo^Q^rms ^ 10 3 cm~ 3 , 
with a peak-to-minimum ratio > 2 for Qo,49 n rms ^ 10 4 cm -3 . As Qonrms is increased, the ionized 
gas becomes concentrated in a thin, dense shell, the peak intensity near the edge rises, and the 
central emission measure changes from EM(0) = 2n^ ms R for small Qo n rms to EM(0) — > (2/3)n^ ms R 

as Qonrms — > CO • 

When t^o S> 1, the present models have the ionized gas concentrated in a thin, dense, shell. 
Figure [5^ shows the ratio of the rms density n rms to the mean density (n) as a function of t^o- The 
highest ionized density occurs at the outer edge of the ionized zone, and Figure [5b shows the ratio 
n e dge/^rms as a function of 7^0 • 

In the low density limit Qo n rms — > 0, the dust optical depth from center to edge t(R) t^q- 
The actual dust optical depth from center to edge is shown in Figure [5b. When the H II region 
develops a dense shell, which occurs for 7^0 > 3, the actual dust optical depth t(R) is significantly 
smaller than the value Tdo- Figure [5b shows that for 7^0 = 40, for example, the actual dust optical 
depth t(R) is only in the range 1-2.3, depending on the values of /3 and 7. 

The shell-like structure is also apparent in the ratio of the peak intensity to the central intensity. 
As seen in Figures H^b-d), dust causes the peak intensity to be off-center. For fixed f3 and 7, the ratio 
of peak intensity to central intensity I(peak)/I(center) increases monotonically with increasing 7^0, 
as shown in Fig. [5](d) . 

Because the shell is dense, radiative recombination is rapid, the neutral hydrogen fraction is 
enhanced, and H atoms can compete with dust to absorb hv > 13.6 eV photons. Figure [6] shows 
/i on , the fraction of the hv > 13.6 eV photons emitted by the star that photoionize H (i.e., are 
not absorbed by dust), as a function of the parameter 7^0 • Results are shown for (3 = 2,3,5 and 
7 = 5, 10, 20. For 2</3<5, 5<7< 20, and < t^ < 40, the numerical results in Figure Eban 
be approximated by the fitting formula 

A»(0,7,t*) - 1 + {2 /3 1 +AB)r d0 + rrB^- (21) 

A - - (99) 

1 + 0.75 7 - 65 /3-°- 44 y ' 

B = l + 0.1( 7 //3) 1 - 5 (23) 
where /3, 7, and t^q are given by ([6]), ([7]), and (fl~6[) . The form of eq. (|21ti23[) has no physical 
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Fig. 5. — For dusty H II regions with 7 = 5, 10, 20 and ft = 2, 3, 5, as a function of Tao" (a) ratio n rm s/ (n) of rms 
density to mean density; (b) ratio n(R)/n lnl s of the edge density to the rms density; (c) center-to-edge dust optical 
depth t(R); (d) ratio of peak emission measure/central emission measure. Results are for {hu)i = 18 eV. 



significance, but eq. (|21|) can be used to estimate the total H ionization rate fionQo in dusty H II 
regions. 

Even for large values of r^o, Fig. [6] shows that ~l/3 of the hy > 13.6 eV pho t ons a re absorbed 
by hydrogen. This contrasts with the uniform-density models of lPetrosian et al.l (|1972l ). where the 
fraction of the hu > 13.6 eV photons that are absorbed by the gas goes to zero as r^ becomes 
large. 
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Fig. 6. — Fraction /j on of the hv > 13.6 eV photons that photoionize H in dusty H II regions with radiation pressure, 
as a function of r^o, for /3 = 2, 3, 5 and 7 = 5, 10, 20. Calculated assuming (hu)i = 18 eV. Th e dotted lines show th e 
fitting formula (|21|l for each case. Also shown is /i Dn calculated for assumed uniform density ( Petrosian et al. 19721 ). 



4. Dust Drift 

4.1. Gas Drag vs. Radiation Pressure 

Eq. (HJ assumes the dust to be tightly coupled to the gas, so that the radiation pressure force 
on the dust can be considered to act directly on the gas. In reality, radiation pressure will drive the 
dust grains through the plasma. If the grains approach their terminal velocities (i.e., acceleration 
can be neglected) then, as before, it can be assumed that the radiation pressure force is effectively 
applied to the gas. However, the motion of the dust grains will lead to changes in the dust/gas ratio, 
due to movement of the grains from one zone to another, as well as because of grain destruction. 
Here we estimate the drift velocities of grains. 

Let Qprtra 2 be the radiation pressure cross section for a grain of radius a. Figure [7] shows 
Q pr (a,A) averaged over blackbody radiation fields with T = 25000 K, 32000 K, and 40000 K, for 
carbonaceous grains and amorphous silicate grains. For spectra characteristic of O stars, (Q pT ) ~ 
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1.5 for 0.02/mi < a < 0.25^m. 




10- 3 0.01 0.1 1 
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Fig. 7. — Spectrum-averaged radiation pressure efficiency factor (Q pr ) as a function of radius, for T = 25000 K, 
32000 K, and 40000 K blackbody spectra. For temperatures characteristic of O stars, (Q pr ) ~ 1-5 to within 20% for 
0.02pm < a < 0.25pm. 



If the magnetic field B = 0, the terminal velocity for a grain at distance r is determined by 
balancing the forces due to radiation pressure and gas drag: 



L(r) 
4?rr 2 c 



ira 2 (Q pi ) 



27ra 2 nkTG(s) 



Vd 



\j2kT jrriYi 



(24) 



where the drag function G(s\ inclu ding both collisional drag and Coulomb drag, can be approxi- 
mated by (jDraine fc Salpeterlll979l ) 



G(s) 
A 



8.s 



, 9tt 2 \ 1/2 feU . , 4 
1 + -tttS 2 + — In A 



30F V" ' 64 

3 kT ( jfeTN 1/2 
2ae |e!7| \7rn e 



(30F/4 + s 3 



6.6 x 10° a 



6„-l 



kT 
~\eU~\ 



1/2 -1/2 



(25) 
(26) 



where £7 is the grain potential, a_5 = a/10 5 cm, and = n e /10 3 cm 3 . The drag force from the 
electrons is smaller than that from the ions by at least -J m e /m p , and can be neglected. The charge 
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on the grains will be determined by coll isional chargi ng and photoelectric emission. Collisional 



charging would result in eU jkT ~ —2.51 (|Spitzerlll968l ). or U ~ — 2.I6T4V. Photoelectric charging 



will dominate close to the star, but is expected to result in potentials U < 10V. Taking \eU /kT\ 
2.5 and In A ~ 14.8 as representative, 



G(s) 



9vr 2 \ 1/2 69.5 
1.50 1 + — s 2 + 



64 J 1 + 453/3^ 



(27) 



Note that G(s) is not monotonic: as s increases from 0, G(s) reaches a peak value ~ 42 for s ~ 0.89, 
but then begins to decline with increasing s as the Coulomb drag contribution falls. At sufficiently 
large s, the direct collisional term becomes large enough that G(s) rises above ~ 42 and continues 
to rise thereafter. 

The drag time for a grain of density p = 3 g cm -3 in H II gas is 

Mv I a_c; \ s , , 

^^ 295 b#j^ yr • (28> 

For n.3 > 0.01 this is sufficiently short that each grain can be assumed to be moving at its terminal 
velocity Vd, with isothermal Mach number s = Vd/\/2kT/mn determined by the dimensionless 
equation 

G(s) = U(y) + Pe-M ^(Q PT ) , y = f . (29) 
L J y z A 

Eq. (129p is solved to find s(r). For 20 < G < 42, there are three values of s for which the drag force 
balances the radiation pressure force. The intermediate solution is unstable; we choose the smaller 
solution] 3 ! which means that s undergoes a discontinuous jump from ~ 0.9 to 6.2 at G ~ 42. The 
resulting terminal velocity v(r) is shown in Figure [5] for 7 values of (Jo,49 n rms- The velocities in the 
interior can be very large, but the velocities where most of the dust is located [r(r) /t(R) > 0.5] 
are much smaller. 

In the outer part of the bubble, where most of the gas and dust are located, the drift velocities 
are much more modest. This is visible in Fig. [8^, where the drift speeds become small asr — > R, but 
is more clearly seen in Fig. [8b, showing drift speeds as a function of normalized optical depth. The 
range 0.5 < T{r)/r{R) < 1 contains more than 50% of the dust, and throughout this zone the drift 
speeds are < 0.3 km s _1 even for Qo,49"rms as large as 10 7 cm~ 3 . With drift speeds Vd < 0.3 km s" 1 , 
grains will not be destroyed, except perhaps by shattering in occasional collisions between grains 
with different drift speeds. However, for large values of n rms , these grains are located close to the 
boundary, the drift times may be short, and the grains may be driven out of the H II and into the 
surrounding shell of dense neutral gas. This will be discussed further below. 



This solution is physically relevant if the drift speed began with s < 0.9 and increased with time. 
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Fig. 8. — Radial drift velocities Vd, r for six different H II regions, all with /3 = 3 and 7 = 10, for T4 = 0.94, 
(hu)i = 18 eV, Q pr = 1.5, \eU/kT\ — 2.5, and B = 0. (a) «d jr vs r/_R. All solutions have large drift velocities near 
the center, which will result in removal of the dust from the interior. Drift velocities increase with increasing Qorirms- 
(b) Vd,r as a function of dust column density r(r). Even if B = or B || r, drift velocities Vd > 75 km s _1 occur only 
in a region with a small fraction of the dust. In most of the volume, drift will not result in grain destruction. 



4.2. Magnetic Fields 

Let e_e = B 2 /167rnkT be the ratio of magnetic pressure to gas pressure. The importance 
of magnetic fields for the grain dynamics is determined by the dimensionless ratio wr drag , where 
to = QB/Mc is the gyrofrequency for a grain with charge Q and mass M in a magnetic field B: 

If \eU/kT\ » 2.5 and In A « 15, then (G(s)/s) » 71 for s < 0.5. 

Let the local magnetic field be B = B(rcos9 + y sin0). The steady-state drift velocity is 



/ F rs , d T drSLg \ / 1 + (wr drag ) 2 cos 2 9 

Vd = ( — ^ — ; V i+(- drag ) 2 ' (31) 

, ^radTdrag\ 1 + (W drag ) 2 COS 2 9 

y d,r = u — — ; , 7— n 2 — . (32) 



M J 1 + (wr dr , 



/ -^radTdrag \ (wT drag ) 2 Sm (? COS 9 

V ^ = {-^) 1 + (- drag ) 2 ' ^ 
/ -^radTdrag \ WT drag sin 

«*. = -(^5-J 1 + (MT ^). • < 34 > 

where v^ r , i^, t^z are the radial and transverse components. If sine? — >■ 0, the magnetic field does 
not affect the radiation-pressure-driven drift velocity, but in the limit sin# — > 1 magnetic effects 
can strongly suppress the radial drift if cjr drag » 1 and cos 9 <C 1. 
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Fig. 9. — Drift timescale tdrift/Qo,49 (see eq.[35j for — 3 and 7 = 5, 10, and 20 (assuming T 4 = 0.94, (hv)i = 18 eV). 
The dust grains are assumed to have (Q pT ) = 1.5 and \eU/kT\ = 2.5. Solid lines are for B = (or B || r). Broken 
lines are for a = 0.03^im, and es<3o,49 = 0.1 and 10 2 . 



The magnetic field strength is uncertain, but it is unlikely that the magnetic energy density 
will be comparable to the gas pressure; hence eb < 0.1. From eq. ([30]) it is then apparent that if 
the magnetic field is strong (ee ~ 0.1), magnetic effects on the grain dynamics can be important 
in low density H II regions, but will not be important for very high densities: (wrd ra g) 2 ^$ 1 for 
n 3 > \7QaZle B - 



4.3. Drift Timescale 

When radiation pressure effects are important, the gas and dust are concentrated in a shell 
that becomes increasingly thin as Qon rras is increased. The drift velocities where most of the dust 
is located are not large (see Fig. Eb), but the grains are also not far from the ionization front. The 
timescale on which dust drift would be important can be estimated by calculating the drift velocity 
at the radius ro.5 defined by r(ro.s) = 0.5t(R). More than 50% of the dust has ro.5 < r < R. 
Figure [9] shows the characteristic drift time 

_ R - r , 5 

Mrift = 7 r • [oO) 

Vd,r{ro.5) 

If no magnetic field is present, the drift velocity depends only on T and the dimensionless quantities 
{cp,T,u,y, (Qpr)} (see eq. |29|) . It is easy to see that for fixed T and Qo^rms; the radius R oc Qq, 
thus tdrift oc Qq. Figure [9] shows £drift/Qo49- For Qo,49 = 1, H II regions with ra rms > 10 3 cm -3 
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have tdrift < 10 6 yr if magnetic effects are negligible. If a magnetic field is present with B 1 r and 
6b = 0.1, then the grain drift is slowed, but drift times of < lMyr are found for n rms > 10 4 cm~ 3 . 
Therefore, compact and ultracompact H II regions around single O stars are able to lower the 
dust /gas ratio by means of radial drift of the dust on time scales of < Myr. However, if the O star 
is moving relative to the gas cloud with a velocity of more than a few km s _1 , then individual fluid 
elements pass through the ionized zone on timescales that may be shorter than the drift timescale, 
precluding substantial changes in the dust-to-gas ratio. 

Grain removal by drift can also occur for giant H II regions. As an example, consider a giant 
H II region ionized by a compact cluster of ~10 3 O stars emitting ionizing photons at a rate 
Qo = 10 52 s _1 . For n rms = 10 3 cm -3 , we have Qo,49 n rms = 10 6 cm -3 , and we see that if B = 0, the 
drift timescale is only thrift ~ 2 x 10 5 yr. If a magnetic field is present with B 1 r and e# = 0.1, 
then from Figure [9] the drift timescale tdrift is increased, but only to ~ 10 6 yr. It therefore appears 
possible for radiation-pressure driven drift to remove dust from giant H II regions provided they 
are sufficiently dense. 

Aside from magnetic effects, the drift speeds at a given location depend only on {Q pr ) and 
T4 (see eq. I29|h Figure [7] shows that (Q pr ) is constant to within a factor ~1.5 for a > 0.010/mi. 
Hence radiation-pressure-drive drift would act to drive grains with a > 0.01/im outwards. Smaller 
grains will drift as well, but more slowly. Because of this, the gas-to-dust ratio in the centers of 
H II regions should in general be lower than the gas-to-dust ratio in the gas prior to ionization. 
The dust-to-gas ratio will first be reduced in the center, where the drift speeds (see Fig. [8]) are 
large. Dust drift will also alter the dust-to-gas ratio in the outer ionized material, initially raising 
it by moving dust outwards from the center. In an initial ly uniform n eutral cloud, the ionization 



front expands rapidly at early times (see, e.g., Fig. 37.3 in lDraindl201ll ) but in gas with 77,3 > 1, at 
late times the ionization front will slow to velocities small enough for dust grains to actually drift 
outward across the ionization front, lowering the overall dust-to-gas ratio within the H II region. 



4.4. Grain Destruction 



Arthur et al.l (120041 ) computed models of uniform density H II regions including the effects of 
dust destruction by sublimation or evaporation, finding that the dust/gas ratio can be substantially 
reduced near the star. If the maximum temperature at which a grain can survive is T su b, and the 
Planck-averaged absorption efficiencies are Q uv and Q u for T = T* and T = T max , then grains will 
be destroyed within a distance r su b with 

% — « 3 (CB 2 (^f(^f <»> 

For parameters of interest (e.g., £39/(30,49 ~ 1, ^39 % 10 2 ) we find r su b/i? s o 1 for n rms < 
10 5 cm -3 , and sublimation would therefore destroy only a small fraction of the dust. 
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As we have seen, we expect radiation pressure to drive grains through the gas, with velocity 
given by eq. (|3ip . Drift velocities > 75km s -1 will le ad to sputteri ng by impacting He ions, with 



sputtering yield Y"(He) ~ 0.2 for 80 < v < 500 km s 1 ( Draindll995l ). For hypersonic motion, the 



grain of initial radius a will be destroyed after traversing a column density 

n HA r =™-^2xlO-a- 5 (IM) cm -s (37) 
n He Y[Re)n \ 0.2 J 

for a grain density p/fi = 1 x 10 23 cm -3 , appropriate for either silicates (e.g., FeMgSi04, ~ 
3.8gcm~ 3 /25mH = 9 x 10 22 cm~ 3 ) or carbonaceous material (2gcm~ 3 /12mn = 1-0 x 10 23 cm -3 ). 
Therefore the dust grain must traverse material with (initial) dust optical depth 

Ar d = a d n R Ar = 0.2a d ^ 2 ici- 5 (^y^) ( 38 ) 

if it is to be substantially eroded by sputtering. However, Fig. 03 shows that even in the absence 
of magnetic effects, v d > 75 km s" 1 occurs only in a central region with < 0.05. Therefore 
sputtering arising from radiation-pressure-driven drift will not appreciably affect the dust content. 



5. Discussion 



5.1. Absorption of Ionizing Photons by Dust 



For a sample of 13 Galactic H II regions, llnoud (|2002l ) used infrared and radio continuum 
observations to obtain the values of /; on shown in Figure [TUJ The estimated values of /i on are much 
larger than would be expected for unifor m H II regions with dust-to-gas ratios comparable to the 
values found in neutral clouds, llnoud (120021 ) concluded that the central regions of these H II regions 
must be dust-free, noting that this was likely to be due to the combined effects of stellar winds 
and radiation pressure on dust. As seen in Fig. [TUl the values of /i on found by Inoue are entirely 
consistent with what is expected for static H II regions with radiation pressure for 5 < 7 < 20 



(corresponding to 0.5 < o d< 



-21 ^ 2), with no need to appeal to stellar winds or grain destruction. 



5.2. The Density-Size Correlation for H II Regions 



H II regions come in many sizes, ranging from H II regions powered by a single O star, to 
giant H II regions ionized by a cluster of massive stars. The physical size of the H II region is 
obviously determined both by the total ionizing output Qq provided by the ionizing stars, and the 
r.m.s. density n rms of the ionized gas, which is regulated by the pressure p e dge of the confining 
medium. With the balance between photoionization and recombination determining the size of an 
H II region, an anticorrelation between size D and densit y n rms is expected, and was observed a s 
soon as large samples of H II regions became available (e.g.. lHabing &: Israellll979l ; lKennicuttlll984l ). 
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Fig. 10. — Photoionizing fraction /i on for 12 Galactic H II regions, as estimated by Inouel 1 20021 ) from infrared and 
radio observations, vs. Qo^gWeTf ' 83 (see text). /i on cannot exceed 1, therefore the high value found for G298.22-0.34 
give some indication of the uncertainties in estimation of /i on . Solid lines: /j on for H II regions with radiation pressure 
for dust characterized by 7 = 5, 10, and 20. Broken line: /i on for uniform H II regions with ad = 10 -21 cm 2 H _1 . 



For dustless H II regions, one expects n rms oc D 1,5 for fixed Qp, but for various samples relations 



close to n rrn s oc D 



2001 



Garav &; Lizano 


1999 


■ 


Kim & Koo 


'gions, 


Kim & Koo 


(2001) attribute 



the n e oc D~ l trend to a "champagne f low" and t he hierarchica l struc ture of the dense gas in the 
star-forming region, but lArthur et al.1 ( 2004 ) and iDopita et al.1 ( 20061 ) argue that the n e oc D^ 1 
trend is a result of both absorption by dust and radiation pressure acting on dust in static H II 
regions. 



Hunt &: Hirashital ([2009) recently reexamined the size-density relationship. They interpreted 
the size-density relation for different observational samples in terms of models with different star 
formation rates [and hence different time evolution of the ionizing output Qo(t)], and differences in 
the density of the neutral cloud into which the H II region expands. Their models did not include 
the effects of radiation pressure on dust; at any time the ionized gas in an H II region was taken to 
have uniform density, resulting in overestimation of the dust absorption. 

Figure [TTh shows a grid of n rms vs. D for the present models, for four combinations of (/3,7). 
While differences between the models with different (/3, 7) can be seen, especially for high Qq and 
high Pedge> the overall trends are only weakly dependent on (3 and 7, at least for 1 < 7 < 20. 

Figure II lb shows the model grid for (3 = 3 and 7 = 5 together with observed values of D and 
n rms from a number of different studies. It appears that observed H II regions - ranging from H II 
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Fig. 11. — Density n rms vs. diameter D. (a) Models with (/3,-y) = (2,1), (2,5), (3,10), and (5,20). Results shown 
were calculated for T4 = 0.94, (hv)i = 18 eV. Solid lines show models with p e d ge fixed, and Qo varying from 10 48 s _1 
to 10 54 s _1 . Broken lines show models with Qo fixed, and p e d g e/fc varying from 10 4 cm -3 K to 10 11 cm -3 K. (b) Model 
grid for /3 = 3, 7 = 10 together wi th observed value s are shown for va rious samples of Ga l actic and extragalactic 
HII regions. Cyan op en triangles: Kennicuttl ( 1984). Bl ue diamonds: Churchwell fc Goss Jl999 | ). G r een c rosses: 
Garav fc Lizano 1 1999h . Cyan cro sses: iKim fc KooT (l200lh. Black open stars: iMartm-Hernandez et al.l J2005h . Red 
solid t riangles: radio sample from Hunt fc Hirashita (120091 ). Red open circles: HST sample from 
(200! 



9). 



Hunt fc Hirashita 



regions ionized by one or at most a few O stars (Qo < 10 50 s 1 ) to "super star clusters" powered 

can be accomodated by the present static 



10 



r.2 



10 



54 , 



by up to 10 3 - 10 5 O stars (Q 
equilibrium models with external pressures in the range 10 4 < p/k < 10 10,3 cm _3 K. Note that 
for diameters D > 10 2 pc, the assumption of static equilibrium is unlikely to be justified, because 
the sound-crossing time (L>/2)/15km s _1 > 3Myr becomes longer than the lifetimes of high-mass 
stars. 



The fact that some H II region samples (e.g., iGarav et al.lll993l ; IKim fc Kod l200ll ) seem to 
obey a n rms oc D^ 1 relationship appears to be an artifact of the sample selection. We see in Fig. 
Illb that the overall sample of H II regions does not have a single n rms -vs.-D relationship. But the 
observations appear to be generally consistent with the current models of dusty H II regions. 
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5.3. Cavities in H II Regions: N49 



Even without dust present, radiation pressure from photoelectric absorption by H and He 
can alter the density profile in a static H II region, lowering the central density and enhancing 
the density near the edge of the ionized region (see Fig. [2h). As seen in Figure H^i, for large 
values of Qo n rms the surface brightness profile can be noticeably flattened. If dust is assumed to be 
present, with properties typical of the dust in diffuse clouds, the equilibrium density profile changes 
dramatically, with a central cavity surrounded by a high-pressure shell of ionized gas pushed out 
by radiation pressure. In real H II regions, fast stellar winds will also act to inflate a low-density 
cavity, or "bubble", near the star; the observed density profile will be the combined result of the 
stellar wind bubble and the effects of radiation pressure. 



The GLIMPSE survey (jChurchwell et al.l I2009TI has disco vered and catalogued numerous in- 
terstellar "bubbles". An example is N49 ([Watson et al.ll2008l ). with a ring of free- free continuum 
emission at 20 cm, surrounded by a ring of 8/im PAH emission. An 06.5V star is located near the 
center of the N49 ring. The image is nearly circularly symmetric, with only a modest asymmetry 
that could be due to motion of the star relative to the gas. The 20 cm image has a ring-peak-to- 
center intensity ratio J(peak)/J(center) ~ 2. 

Is the density profile in N49 consistent with wh at is expected for ra diation pressure acting on 
dust? From the 2 89 Jy flux from N49 at A = 20 cm (jHelfand et alJbood ) and distance 5.7±0.6kpc 
(jChurchwell et al.1 120061 ). the stellar source has Qo,49 ~ (0.78 ± 0.16) / /ion- If /ion ~ 0.6, then 
Qo,49 ~ (1-3 ± 0.3). The H II region, with radius (0.018 ± 0.02) deg, has n rms » 197 ± 63cm~ 3 . 
Hence Qo,49"<rms ~ 260 cm~ 3 . If <7d-2i = 1 ; then Tdo ~ 1-3. From Fig. [6^ we confirm that /; on « 0.6 
for T d0 « 1.3. 

Figure [5ti shows that an H II region with t^o = 1-3 is expected to have a central minimum in 
the emission measure, but with J(peak)/ /(center ) ~ 1.3 for /3 = 3, 7 = 10, whereas the observed 
J(peak)//(center) ~ 2. The central cavity in N49 is therefore significantly larger than would be 
expected based on radiation pressure alone. While the effects of radiation pressure are not negligible 
in N49, the observed cavity must be the result of the combined effects of radiation pressure and a 
dynamically-important stellar wind (which is of course not unexpected for an 06.5V star). 



5.4. Lyman-o 

The original ionizing photon deposits a radial momentum hvi/c at the point where it is ab- 
sorbed by either a neutral atom or a dust grain. A fraction (1 — /km) of the ionizing photons are 
absorbed by dust; this energy is reradiated isotropically, with no additional force exerted on the 
emitting material. Because the infrared optical depth within the H II region is small, the infrared 
emission escapes freely, with no dynamical effect within the H II region. 

A fraction /j on of the ionizing energy is absorbed by the gas. Subsequent radiative recombi- 
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nation and radiative cooling converts this energy to photons, but the isotropic emission process 
itself involves no net momentum transfer to the gas. We have seen above that the H II can have 
a center-to-edge dust optical depth t(R) ~ 1.6 for t^o > 5, or Qo,49"rms ^ 10 2 cm -3 (cf. Fig. 
[5^ with /3 = 3, 7 = 10). This optical depth applies to the hv > 13.6 eV ionizing radiation; the 
center-to-edge optical depth for the hv < 3.4 eV Balmer lines and collisionally-excited cooling lines 
emitted by the ionized gas will be significantly smaller, and much of this radiation will escape dust 
absorption or scattering within the H II region. That which is absorbed or scattered will exert a 
force on the dust at that point only to the extent that the diffuse radiation field is anisotropic. We 
conclude that momentum deposition from the Balmer lines and collisionally-excited cooling lines 
within the ionized zone will be small compared to the momentum deposited by stellar photons. 

Lyman-a is a special case. At low densities (n <C 10 3 cm -3 ), ~ 70% of Case B recombinations 
result in emission of a Ly-a pho ton, increasing to > 95% f or n > 10 5 cm -3 as a result of collisionally- 
induced 2s — > 2p transitions (jBrown &: Mathews Il970h . After being emitted isotropically, the 
photon may scatter many times before either escaping from the H II or being absorbed by dust. 
Most of the scatterings take place near the point of emission, while the photon frequency is still 
close to line-center. On average, the net radial momentum transfer per emitted photon will likely 
be dominated by the last scattering event before the photon escapes from the H II region, or by the 
dust absorption event if it does not. At a given point in the nebula, the incident photons involved 
in these final events will be only moderately anisotropic. Since there is less than one Ly-a photon 
created per case B recombination, the total radial momentum deposited b y these final events will be 
a small fraction of the radial momentum of the original ionizing photons. Henney &i Arthur (jl998l ) 
estimate that dust limits the Ly-a radiation pressure to ~6% of the gas pressure. We conclude 
that Ly-a has only a minor effect on the density profile within the ionized zone. 



5.5. H II Region Expansion 

H II regions arise when massive stars begin to emit ionizing radiation. The development of the 
H II region over time depends on the growth of the ionizing output fro m the centra l star, and the 



expansion of the initially-high pressure ionizing gas. Many authors (e.g.. lKahnlll954l ; ISpitzerlll978l ) 
have discussed the development of an H II region in gas that is initially neutral and uniform. If the 
ionizing output from the star turns on suddenly, the ionization front is initially "strong R-type", 
propagating supersonically without affecting the density of the gas, slowing until it becomes "R- 
critical" , at which point it makes a transition to "D-type" , with the ionization front now preceded 
by a shock wave producing a dense (expanding) shell of neutral gas bounding the ionized region. 

While the front is R-type, the gas density and pressure are essentially uniform within the 
ionized zone. When the front becomes D-type, a rarefaction wave propagates inward from the 
ionization front, but the gas pressure (if radiation pressure effects are not important) remains 
relatively uniform within the ionized region, because the motions in the ionized gas are subsonic. 
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When radiation pressure effects are included, the instantaneous density profile interior to the 
ionization front is expected to be similar to the profile calculated for the static equilibria studied 
here. Let Vi be the velocity of the ionization front relative to the star. When the ionization front is 
weak D-type, the velocity of th e ionization front relative to the ionized gas just inside the ionization 
front is ~ 0.5V^ (|Spitzerlll978l ). Given the small dust drift velocities v^ T near the ionization front 
(i.e., r(r) — > t(R) in Fig. [8]), dust is unable to drift outward across the ionization front as long as the 
ionization front is propagating outward with a speed (relative to the ionized gas) Vi > 0.1km s _1 



6. Summary 

1. Dusty H II regions in static equilibrium consist of a three-parameter family of similarity 
solutions, parametrized by parameters /?, 7, and a third parameter, which can be taken to be 
Qo,49^rms or T^o (see eq. [To]) . The j3 parameter (eq. [6]) characterizes the relative importance 
of hv < 13.6 eV photons, and 7 (eq. [7]) characterizes the dust opacity. A fourth parameter - 
e.g., the value of n rms or Qo,49 ~ determines the overall size and density of the H II region. 

2. Radiation pressure acting on both gas and dust can strongly affect the structure of H II 
regions. For dust characteristic of the diffuse ISM of the Milky Way, static H II regions with 
Qo,49^rms ^ 10 2 cm~ 3 will have nearly uniform density, but when Qo,49^rms 3> 10 2 cm~ 3 , 
radiation pressure acts to concentrate the gas in a spherical shell. 

3. For given (5 and 7, the importance of radiation pressure is determined mainly by the parameter 
Tdo (see eq. [TBI) . When 7^0 > 1, radiation pressure will produce a noticeable central cavity. 

4. If the dust-to-gas ratio is similar to the value in the Milky Way, then compression of the 
ionized gas into a shell limits the characteristic ionization parameter: U1/2 ^5 0.01, even for 
Qo"rms > 1 (see Fig. [3]). 

5. For Tdo > 1, compression of the gas and dust into an ion ized shell leads to a substantial increase 



(compared to the estimate by iPetrosian et al.l Il972l ) in the fraction /; on of hv > 13.6 eV 



photons that actually ionize H, relative to what would have been estimated for a uniform 
density H II region, as shown in Fig. [6j Eq. (f2Tj) allows f lon to be estimated for given Qo^rms, 
(3, and 7. Galactic H II regions appear to have values of /; on consistent with the present 
results for H II regions with radiation pressure (see Fig. [TOj) . 

6. Interstellar bubbles surrounding O stars are the result of the combined effects of radiation 
pressure and stellar winds. For the N49 bubble, as an example, the observed ring-like free- 
free emission profile is more strongly peaked than would be expected from radiation pressure 
alone, implying that a fast stellar wind must be present to help create the low-density central 
cavity. 

7. For static H II regions, dust drift would be important on time scales ^ llVlyr for (5o,49^Tms ^ 
10 3 cm -3 . Real H II regions are not static, and the dust will not drift out of the ionized gas 
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because the ionization front will generally be propagating (relative to the ionized gas just 
inside the ionization front) faster than the dust drift speed < 1km s _1 (see Fig. [8]). 

I am grateful to Bob Benjamin and Leslie Hunt for helpful discussions, to R.H. Lupton for making 
available the SM graphics package, and to the anonymous referee for suggestions that improved the 
paper. This research made use of NASA's Astrophysics Data System Service, and was supported 
in part by NASA through JPL contract 1329088, and in part by NSF grant AST 1008570. 
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